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Abstract 

Nonsingular,  stressed,  dislocation  (wall)  profiles  are  shown  to  be  1-d  equilibria  of  a  non-equilibrium 
theory  of  Field  Dislocation  Mechanics  (FDM).  It  is  also  shown  that  such  equilibrium  profiles 
corresponding  to  a  given  level  of  load  cannot  generally  serve  as  a  traveling  wave  profile  of  the  governing 
equation  for  other  values  of  nearby  constant  load;  however,  one  case  of  soft  loading  with  a  special  form 
of  the  dislocation  velocity  law  is  demonstrated  to  have  no  ‘Peierls  barrier’  in  this  sense.  The  analysis  is 
facilitated  by  the  formulation  of  a  1-d,  scalar,  time-dependent,  Hamilton-Jacobi  equation  as  an  exact 
special  case  of  the  full  3-d  FDM  theory  accounting  for  non-convex  elastic  energy,  small,  Nye-tensor 
dependent  core  energy,  and  possibly  an  energy  contribution  based  on  incompatible  slip.  Relevant 
nonlinear  stability  questions  are  formulated  in  a  non-equilibrium  setting.  Elementary  averaging  ideas 
show  a  singular  perturbation  structure  in  the  evolution  of  the  (unsymmetric)  macroscopic  plastic 
distortion,  thus  pointing  to  the  possibility  of  predicting  generally  rate-insensitive  slow  response 
constrained  to  a  tensorial  ‘yield’  surface,  while  allowing  fast  excursions  off  it,  even  though  only  simple 
kinetic  assumptions  are  employed  in  the  microscopic  FDM  theory.  The  emergent  small  viscosity  on 
averaging  that  serves  as  the  small  parameter  for  the  perturbation  structure  is  a  robust,  almost-geometric 
consequence  of  large  gradients  of  slip  in  the  dislocation  core  and  the  presence  of  a  large  number  of 
dislocations  in  the  averaging  volume.  In  the  simplest  approximation,  the  macroscopic  yield  criterion 
displays  anisotropy  based  on  the  microscopic  dislocation  line  and  Burgers  vector  distribution,  a 
dependence  on  the  Laplacian  of  the  incompatible  slip  tensor  and  a  nonlocal  term  related  to  a  Stokes- 
Helmholtz-curl  projection  of  an  ‘internal  stress’  derived  from  the  incompatible  slip  energy. 

1.  Introduction 

We  are  concerned  here  with  the  development  of  an  idealized  mathematical  description  of 
nominally  crystalline  solids  containing  dislocations  as  the  only  defects.  The  interest  is  in 
behavior  at  the  time  scale  of  ~  microseconds  or  more  and  therefore  in  theory  that  can  be 
routinely  used  for  attempting  to  understand  such  behavior.  It  is  our  belief  that  an  appropriate, 
nonlocal  li eld/continuum  description  that  involves  fields  that  are  smooth  over  a  length  scale  of 
interatomic  length  and  larger  can  be  adequate  for  such  a  purpose.  Kinetic  energy  of  atomic  and 
subatomic  vibrations  below  this  length  scale,  and  time  averaged  over  periods  of  microseconds,  is 
characterized  as  dissipation  within  such  a  framework.  Thus,  in  concept,  we  exchange  the  atoms 
and  their  (sub)femtosecond  vibration  periods  in  Molecular  Dynamics  with  the  dissipative 
evolution  of  a  mathematical  dislocation  density  field  (capable  of  representing  ‘individual’ 
dislocations)  within  a  nominally  elastic  continuum  characterized  eventually  by  Cauchy-Born 
crystal  elasticity  (Milstein,  1982;  Ortiz  and  Phillips,  1999).  The  stored  energy  of  the  body  is 
augmented  with  some  additional  contributions  to  represent  nonlocal  energy  content  that  cannot 
be  ascribed  to  local  elastic  response,  but  hopefully  with  enough  physical  meaning  so  as  to  be 
capable  of  unambiguous  definition  from  quantum/atomistic  studies.  Our  main  goal  in  pursuing 
this  direction  of  work  is  really  the  prediction  of  the  dynamics  of  mesoscale  microstructure  and  its 
effect  on  macroscopic  properties;  since  temporal  coarse-graining  of  MD  remains  elusive  and  a 
rational  discrete-to-continuum  transition  in  the  3-d,  time-dependent  setting  even  more  so  in  the 
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presence  of  large  numbers  of  defects,  we  believe  that  it  is  important  to  understand  exactly  what 
can  and  cannot  be  achieved  by  an  atomic  scale  field  theory  grounded  in  sound  kinematics, 
conservation  statements,  and  the  possibility  of  fitting  necessary  physics  from  subatomic  scales, 
since  such  a  framework  has  the  potential  of  being  an  effective  unified  tool  for  achieving  the 
desired  scale  transitions.  We  provide  a  very  first  demonstration  of  this  overall  approach  in  the 
paper,  and,  of  course,  much  remains  to  be  done. 

We  pursue  work  here  in  the  geometrically  linear  setting  for  two  reasons.  First,  the  continuum 
mechanics  and  physical  ideas  involved  are  sufficiently  non-standard  to  warrant  leaving  aside  the 
complications  due  to  geometric  nonlinearity  for  the  moment.  Second,  for  a  part  of  the  present 
model,  the  extension  to  finite  defonnations  exists  in  Acharya  (2004),  and  extending  that  earlier 
work  to  account  for  the  new  developments,  while  technically  intricate  from  the  continuum 
mechanics  perspective,  would  be  conceptually  straightforward. 

An  interesting  conceptual  issue  comes  up  in  this  work  related  to  the  definition  of  the  condition 
for  the  onset  of  motion  of  a  single  dislocation,  especially  when  compared  to  the  treatment  of  the 
same  question  as  initiated  in  Peierls  (1940)  based  on  the  static,  classical,  elastic  theory  of 
dislocations.  Due  to  the  translational  invariance  of  the  sum  of  the  linear  elastic  and  the 
interplanar  energy  of  an  infinite  body  in  that  analysis,  the  dislocation  core  could  be  situated 
anywhere  on  the  slip  plane  with  the  same  energy  cost.  This  prompted  Peierls  to  conclude  that, 
within  that  theory,  there  could  not  be  a  threshold  stress  to  initiate  rigid  motion  of  a  dislocation  to 
move  it  from  one  location  to  another.  A  stronger  result  is  shown  by  Movchan  et  al.  (1998)  in  that 
an  equilibrium  solution  to  the  problem  with  a  non-zero,  but  small,  applied  stress  does  not  exist. 
Peierls’s  (1940)  conclusion  and  preliminary  work  prompted  the  accounting  of  lattice  discreteness 
in  the  theory  of  the  Peierls  stress  (Nabarro,  1947)  followed  by  a  rich  literature  (e.g.  Picu,  2002), 
whereby  invariance  of  the  energy  with  respect  to  all  continuous  spatial  translations  of  its 
argument  fields  is  exchanged  for  discrete  invariance  respecting  the  symmetry  of  the  lattice  (e.g. 
Lu  et  al.,  2000). 

The  model  presented  in  this  paper  is  translationally  invariant.  It  is  also  dynamic,  even  in  the 
absence  of  inertia.  In  the  context  of  the  static  approach  to  the  Peierls  stress  problem,  while  the 
absence  of  an  equilibrium  configuration  under  small  load  is  a  compelling  result,  it  is  perhaps  fair 
to  say  that  this  result  by  itself  does  not  say  anything  about  how  the  dislocation  profile  evolves 
under  non-zero  load  -  e.g.  does  instability  lead  to  a  complete  collapse  of  the  slip  configuration, 
or  a  very  severe  distortion  with  no  semblance  to  a  slip  profile  representing  a  dislocation,  or  does 
it  mean  rigid  translation  of  the  equilibrium  profile  corresponding  to  no  load?  Or,  very 
importantly,  does  the  dislocation  density  profile  evolve,  but  with  indiscernible  translation  of  the 
dislocation  core  up  to  a  certain  value  of  the  load  after  which  there  is  an  almost  rigid  translation? 
What  effect  does  rate  of  loading  and  temperature  have  on  this  experiment?  It  seems  that  a 
dynamical  model,  with  a  basis  in  mechanics  (e.g.  the  discrete  Peierls  model  of  Movchan  et  al., 
1998),  is  required  for  such  a  purpose.  Thus,  we  pose  the  question  of  existence  of  a  barrier  to 
rigid  translation  as  one  of  determining  when  dislocation-like  equilibrium  profiles  under  zero  load 
may  not  serve  as  traveling  wave  profiles  under  load,  realizing  full  well  that  the  practically 
relevant  criterion  would  be  related  to  “almost  rigid’  translation,  but  defining  the  latter  precisely 
is  far  more  delicate  than  the  former.  In  a  simplified  1-d  setting,  we  show  that  the  answer  varies 
with  the  choice  of  dislocation  velocity  law  and  loading  mode  (displacement  or  traction  control), 
and  clearly  translational  invariance  is  not  the  crux  of  the  issue  in  this  model.  Of  course,  our 
general  theory  does  not  preclude  the  inclusion  of  discrete  translational  invariance  in  the  model 
through  physically  appropriate  modifications  to  the  form  of  the  energy  density  function,  in  the 
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spirit  of  the  semi-discrete  Peierls  model  (Lu  et  al.,  2000)  and/or  through  a  restatement  in  the 
language  of  discrete  calculus  (e.g.  Ariza  and  Ortiz,  2005). 

This  paper  is  organized  as  follows:  in  Section  2  we  settle  on  some  notational  conventions. 
Section  3  describes  the  theoretical  framework.  In  Section  4  a  simple,  but  exact,  model  problem  is 
developed  and  the  general  theory  is  applied  to  it.  We  obtain  some  special  exact  solutions  for 
equilibria  and  a  travelling  wave,  and  formulate  some  questions  of  dynamical  stability  whose 
answers  would  be  of  direct  physical  interest.  In  Section  5  we  consider  preliminary  implications 
of  our  atomic-scale  model  on  meso  and  macroscopic  response.  In  Section  6  we  present  some 
comparative  observations  on  aspects  of  our  model  with  time-dependent  Ginzburg  Landau  models 
and  level-set  propagation. 

2.  Notation 

A  superposed  dot  on  a  symbol  represents  a  material  time  derivative.  The  statement  a b  is 
meant  to  indicate  that  a  is  being  defined  to  be  equal  to  b  .  The  summation  convention  is  implied. 
We  denote  by  Ab  the  action  of  the  second-order  (third-order,  fourth-order)  tensor  A  on  the 
vector  (second-order  tensor,  second-order  tensor)  b ,  producing  a  vector  (vector,  second-order 
tensor).  A  •  represents  the  inner  product  of  two  vectors,  a  :  represents  the  trace  inner  product  of 
two  second-order  tensors  (in  rectangular  Cartesian  components,  A.B  =  A.  Btj )  and  matrices.  The 

symbol  AB  represents  tensor  multiplication  of  the  second-order  tensors  A  and  B  .  The 
curl  operation  and  the  cross  product  of  a  second-order  tensor  and  a  vector  are  defined  in  analogy 
with  the  vectorial  case  and  the  divergence  of  a  second-order  tensor:  for  a  second-order  tensor  A  , 
a  vector  v  ,  and  a  spatially  constant  vector  field  c , 

(^xv)rc  =  ^Arcjxv  Vc 

(divAf  c  =  div(^ATcj  Vc  (1) 

(curl  Af  c  =  curl  i^AJ cj  Vc. 

In  rectangular  Cartesian  components, 

(Axv).m=emjkAyVk 

(divA)i  =  Ayj  (2) 

iCUrlA)im  =emjAk,p 

where  emjk  is  a  component  of  the  third-order  alternating  tensor  X .  The  spatial  derivative,  for 

the  component  representation  is  with  respect  to  rectangular  Cartesian  coordinates  on  the  body. 
For  all  manipulations  with  components,  we  shall  always  use  such  rectangular  Cartesian 
coordinates.  The  symbol  div  represents  the  divergence,  grad  the  gradient,  and  div  grad  the 

Laplacian.  We  often  have  occasion  to  use  the  identity  curl  curl  (•)  =  grad  div  (•)- div  grad  (•) , 

often  for  an  argument  for  which  div(-  )  =  0 

3.  Field  Dislocation  Mechanics 

The  theory  utilizes  a  continuum  description  of  defect  density  based  on  the  incompatibility  of 
elastic  distortion.  As  such,  it  is  potentially  applicable  to  materials  for  which  an  elastic  response  is 
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defined  (e.g.  amorphous  materials  like  metallic  glasses);  in  this  paper  we  associate  the  defect 
density  field  a  with  a  continuously  distributed  dislocation  density  in  crystalline  materials  with 
single  dislocations  viewed  as  non-singular,  spatially  localized  structures  in  this  continuously 
distributed  density  field.  The  utility  of  such  a  representation,  even  for  the  visualization  of  core 
structures  from  computer  simulations,  is  demonstrated  in  Hartley  and  Mishin  (2005)1. 

The  elastic  distortion  of  the  lattice  is  represented  by  the  (generally  nonsymmetric)  tensor  Ue . 
Due  to  lattice  incompatibility,  the  elastic  distortion  is  not  generally  a  gradient  of  a  vector  field 
and  we  write  its  incompatibility  field  as 

curlUe—.a.  (3) 

Since  we  insist  that  the  total  observable  displacement  gradient  field  (from  an  arbitrarily  chosen 
reference  configuration,  say  B)  be  the  sum  of  the  plastic  distortion,  Up  (with  respect  to  that 
reference),  and  the  elastic  distortion,  the  plastic  distortion  is  also  not  in  general  a  gradient; 

grad  u  =  Ue  +  UP  =>  curl  Ue  =  —curl  Up .  (4) 

It  is  written  as  a  sum  of  a  gradient  and  an  incompatible  part  (that  cannot  be  expressed  as  a 
gradient) 

Up  =  grad  z-%  ■  (5) 

The  incompatible  part  results  from  the  distribution  a  through  the  fundamental  geometrical 
equation  of  incompatibility 


with  the  side  condition 


to  ensure  that  when  a 


curl Up  =  -a  =>  curl x~a  (6) 

div  x  =  0  on  B 

(7) 

Xn  =  0  on  the  boundary  SB  with  unit  normal  n 
=  0  the  incompatible  part  x  vanishes  identically  on  the  body.  Thus,  (5) 


may  be  interpreted  as  the  Stokes-Helmholtz  decomposition  of  Up  into  a  gradient  field  and  a  curl 
field  (see  (22)  below).  The  compatible  part  grad  z  depends  upon  the  history  of  plastic  straining 
and  accumulates  the  compatible  increments  of  the  plastic  strain  produced  by  the  motion  of  the 
dislocation  density  through  the  equation 

div  grad  z  =  div  ( a  x  V ) .  (8) 

In  this  model  of  dislocation  mechanics,  the  total  displacement  field,  u ,  does  not  represent  the 
actual  physical  motion  of  atoms  involving  topological  changes  but  only  a  consistent  shape 
change  and  hence  is  not  required  to  be  discontinuous.  However,  the  stress  produced  by  these 
topological  changes  in  the  lattice  is  adequately  reflected  in  the  theory  through  the  utilization  of 
incompatible  elastic/plastic  distortions.  As  usual  in  continuum  plasticity,  the  stress  is  a  function 
of  the  elastic  distortion  (in  the  linear  elastic  case  given  by  T  =  CUe).  While  not  necessary  in  the 
context  of  the  general  theory,  in  this  paper  we  focus  on  defect  dynamics  in  the  presence  of  static 
force  balance  characterized  by 

divT  =  0 .  (9) 

The  definition  (3)  locally  renders  a  as  a  density  of  lines  carrying  a  vectorial  attribute  (the 
Burgers  vector),  and  it  is  natural  to  associate  a  velocity  field  that  moves  these  line  segments  in 


1  The  paper  contains  some  unfortunate  typographical  errors  that  cloud  the  exposition  of  the  method.  Also,  the  link 
that  the  authors  wish  to  demonstrate  with  Bilby’s  lattice  correspondence  functions  could  not  be  followed  by  this 
author,  but  thankfully  that  link  is  not  a  relevant  part  of  the  essential  construction  of  the  paper. 
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the  body.  This  idea  can  be  made  geometrically  rigorous  and  results  in  the  fundamental 
conservation  law  (due  to  Mura  (1963)  who  attributes  it  to  Kroner): 

a  =  -curl  (axf),  (10) 


where  the  field  V  at  any  spatio-temporal  location  represents  the  velocity  of  the  infinitesimal 
dislocation  segment  at  that  location.  Gathering  all  equations,  the  complete  theory  reads  as 

curl  %  =  a  | 


div  %  =  0 

div  [grad  z)  =  div  (axK) 
Ue  :=  grad  (u  —  z)  +  X 


div 


T  U( 


a  =  — curl  (orxK) 


on  B. 


(11) 


As  for  boundary  conditions, 


Xn  =  o 

[grad  z  —  axV)n  =  0 


on  dB 


(12) 


are  imposed  along  with  standard  conditions  on  displacement  and/or  traction.  For  the  dislocation 
density  field,  analysis  from  the  linear  partial  differential  equation  point  of  view  indicates 
(Acharya,  2003)  that  it  suffices  to  prescribe  a[V  ■«)  on  inflow  parts  of  the  boundary,  but  we 


believe  that  the  nonlinear  problem  admits  other  physically  motivated  possibilities  (Acharya  and 
Roy,  2006). 


3.1  Dissipation 


In  (11 )-( 12)  the  functions  T  and  V  need  to  be  constitutively  specified.  We  seek  guidance  in 
doing  so  by  following  a  global  analog  of  parts  of  a  procedure  initiated  by  Coleman  and  Noll 
(1963)  and  Coleman  and  Gurtin  (1967)  for  local  constitutive  equations.  A  specific  free  energy 
density  function  is  introduced  as  follows: 


Y  =  'v{uesym)  +  s{\a\l)G{x)  +  ^a:cc, 
1 

s  (x)  :=  xm ,  x  >  0,  m  1, 


(13) 


where  /  is  a  typical  interatomic  distance,  e.g.  Burgers  vector  magnitude,  of  the  crystal. 

The  first  term  contains  the  contribution  to  the  stored  energy  due  to  deformation  of  the  lattice 
due  to  the  application  of  loads  and  the  presence  of  dislocations  -  since  we  are  interested  in 
nonlinear  crystal  elasticity  respecting  the  symmetries  of  a  lattice,  this  contribution  is  in 
inherently  nonlinear  due  to  lattice  periodicity. 

The  second  term  reflects  a  local  dependence  on  incompatible  slip  in  dislocated  regions.  To  the 
extent  that  any  energy  is  physically  understood,  it  seems  reasonable  to  demand  that  it  be  defined 
purely  in  terms  of  the  observed  state.  It  is  for  this  reason  that  only  the  incompatible  slip, 
completely  identifiable  from  the  observable,  instantaneous  dislocation  density  field,  occurs  as  an 
argument  of  the  stored  energy  and  not  the  whole  of  the  plastic  distortion  Up  as  the  latter 
contains  the  compatible  part,  grad  z ,  whose  knowledge  requires  knowledge  of  the  history  of 
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dislocation  motion  ( 1 1 )  3 .  The  function  s  ensures  that  the  energy  due  to  incompatible  slip  is 
taken  into  account  only  in  regions  where  the  lattice  has  a  ‘fault’,  as  reckoned  by  the  dislocation 
density.  Thus,  even  though  a  volumetric  density,  if  dislocated  regions  are  localized  around 
planes,  this  energy  contribution  is  also  localized  similarly. 

The  third  term,  with  a  =  O  ( /// 2  j  where  //  is  a  typical  elastic  modulus,  represents  the  fact  that 


the  dislocation  density  in  the  body  cannot  be  greater  in  order  of  magnitude  than  roughly  ~  1// . 
The  physical  justification  for  the  statement  is  that  a  gradient  in  elastic  distortion  can  be  reliably 
measured  over  a  minimal  distance  of  /  in  a  crystal.  Thus,  a  is  to  be  thought  of  as  a  small 
parameter.  Alternatively,  the  term  may  also  be  thought  of  as  an  energy  contribution  in  the 
material  that  arises  in  the  presence  of  high  gradients  of  incompatible  lattice  distortion. 

Taken  together,  the  last  two  terms  in  (13)  are  meant  to  represent  the  fact  that  a  solely  local 
elastic  energy  response,  even  though  nonconvex,  may  not  suffice  to  represent  the  energy  content 
of  a  material  element  at  the  atomic  scale,  if  it  happens  to  contain  incompatible  gradients  of  slip. 
Clearly,  we  are  thinking  of  the  body  as  a  continuum  containing  material  points  even  in  the 

regions  between  successive  atomic  planes  and  that  such  points  are  capable  of  transmitting  forces 

between  themselves.  We  hope  that  the  dominant  qualitative  characteristics  of  actual  nonlinear 
and  nonlocal  energy  of  a  defected  lattice  can  be  represented  by  the  postulated  general  energy  in 
(13)  but,  of  course,  proving  the  correspondence  in  a  physically  rigorous  manner  is  not  obvious, 
and  less  so  mathematically. 

We  define  the  mechanical  dissipation  in  the  body  as  the  difference  of  the  power  of  the  applied 
forces  and  rate  of  change  of  the  stored  energy,  i.e. 

D—  /  ( Tn)-uda —  I  i i/dv  —  /  T  \gradiidv—  /  y/  dv ,  (14) 

Jffl'  '  J  B  J  B  J  B 


where  we  have  ignored  inertia  represented  in  the  displacement  field,  but  including  it  poses  no 
special  problems. 

Our  modest  goal  now  is  to  make  choices  for  V  and  T  that  ensure  D  >  0 .  On  assuming 

T  _ 


dU. 


(15) 


sym 


consistent  with  elasticity,  the  dissipation  may  be  written  as 


D  =  J  T:Updv  —  J  sa:ddv  —  J  %dv  — J  Gr^,/ (|«|/)/ |ar|t/v . 


(16) 


V  needs  to  appear  in  (16)  in  order  to  identify  its  driving  force.  To  this  end,  we  first  show  that  the 
structure  of  the  equations  ( 1 1)-(12)  imply  Up  =  a  x  V  .  Equations  (5),  (1 1 )  , 3  6  imply 


curl  Up  =  curl  (or  x  F) 
div  Up  =  div  (a  x  V) . 


(17) 


(17),  implies  the  existence  of  a  vector  field2,  say  a  ,  such  that  Up  —  axV  =  grad  a  ,  and  then 


(17)  2  implies  that  a  satisfies  Laplace’s  equation 

div  grad  a  =  0  on  B  .  (18) 

Now,  on  the  boundary  3B ,  Upn  =  —%n  +  grad  zn  —  (axV)n  from  (5),(12) , , ;  therefore, 


2  Assuming  the  required  degree  of  triviality  in  the  topology  of  the  body. 
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(, grada)n  =  0  on  d B.  (19) 

Consequently,  (18)  and  (19)  along  with  the  uniqueness  of  solutions  Laplace’s  equation  with 
Neumann  boundary  conditions  (up  to  constant  fields)  implies 

grad  a  —  0 


Up  =axV 


Hence,  the  dissipation  may  be  written,  utilizing  (5),  (1 1)6  and  (20)  as 
D  =  f  T :  [a x  V) dv  +  J  sa :  curl (axV)dv—  I,R  ■%dv  ; 


R  :=  s[\a\i)^- . 


dx 


(20) 


(21) 


In  writing  (21)  we  have  essentially  ignored  the  contribution  from  the  last  term  in  (16).  For  m 
large,  s  is  essentially  the  constant  function  with  value  1,  except  near  zero.  By  carrying  along  the 
ignored  contribution  through  to  equation  (27)  below,  it  can  be  seen  that  ignoring  this  tenn  from 
(21)  onwards  implies  no  loss  in  generality.  In  any  case,  this  whole  exercise  is  meant  to  motivate 
a  constitutive  assumption  for  V ,  and  ignoring  the  said  tenn  may  just  as  well  be  considered  as 
part  of  this  assumption. 

To  manipulate  the  last  term  on  the  right-hand- side  of  (21),  we  utilize  a  Stokes-Helmholtz 
resolution  of  a  square-integrable  tensor  field  with  square-integrable  gradients  due  to  Friedrichs 
(see,  e.g.,  Jiang,  1998,  Theorems  5.8,  5.2).  The  resolution  states  that  given  such  a  tensor  field  R  , 
there  exists  a  unique  (up  to  a  constant)  vector  field  gR  and  a  unique  tensor  field  WR  satisfying 


—div  grad  WR  =  curl  R 
WRxn  —  0  on  dB 
div  grad  gR  =  div  R  ; 
R  =  curl  WR  +  grad  gR 


div  WR  =  0 


on  B 


(22) 


[grad  gR—  R)n  —  0  on  dB  and 
on  B. 

In  keeping  with  the  decomposition  of  the  plastic  distortion,  we  alternatively  refer  to  curl  WR  as 

Xr  :=  curlWR .  (23) 

Utilizing  this  decomposition,  (1 1)  2  6  and  (12)  j , 

J"  R:%dv  =  —J  WR:curl  (axV)dv .  (24) 

Therefore, 


D=/Xr  +  [scurl  a  +  curl  WR  ]r  j  aj  •  V  dv  oc '■[(axV'jxn 


da 


(25) 

utilizing  the  fact  that  WR:\L(axV)xn[  =  —[axV):(WRxn)  =  0  on  dB  due  to  (22). 

We  assume  that  the  boundary  tenn  is  dominated  by  the  interior  term,  thus  identifying  the 
‘ driving  force  ’  for  the  dislocation  velocity  field  as 

V  — >  x[t  +  [scurla  +  Xr\\<*  •  (26) 

Here  we  consider  simple  kinetic  assumptions  of  the  form 


3  In  applying  this  decomposition  to  (5),  (7)  note  that  W  p  xn  =  0  on  dB  implies  that  the  line  integral  of  W  p 
on  arbitraiy  closed  curves  on  the  boundary  dB  vanishes  so  that  by  Stokes  theorem,  curl  n  =  —  %tl  =  0  on 
dB. 
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X  ( T  +  [scurl  a  +  %R  f  )a 
V  =  — ^ ^ - !—  ;  B>  0 

B 

>  [  B 

B  =  \  .  ,  where  B,  B*  are  material  constants, 

[B  \a\ 


(27) 


and  the  physical  dimensions  of  B  is  that  of  stress  x  (length)  x  (velocity)  * .  It  is  to  be  noted 
that  a  ‘classical’  dislocation  is  identified  in  this  model  as  a  region  of  linear  dimension  of  /  (the 
core)  over  which  the  incompatible  slip  varies  continuously  by  an  0(1)  amount;  assumption  (27) 


implies  that  in  this  model  the  dislocation  velocity  field  actually  varies  significantly  in  this  region. 

With  this  constitutive  assumption,  the  governing  equation  for  a ,  a  tensorial  conservation  law, 
is  given  by 


a  =  —curl 


a  x 


X  \T  +  [ scurl  a  +  Zr  f ) a 


H\a\) 


(28) 


The  corresponding  governing  equation  for  the  plastic  distortion  is  inferred  from  (20)  as 

X  [t  +  \—s  div  grad  %  +  %R  f )  curIUp 

Up  =  curlUp  x  — ^ - ~c~Y\ - ix - 1 - •  (29) 

B  \jcurlUp  |) 

Of  course,  the  jump  conditions/continuity  requirements  for  the  plastic  distortion  cannot  be 
inferred  from  (29)  without  knowledge  of  the  form  of  the  physical  conservation  law  (28).  Equation 
(29)  is  a  ‘nonlocal’  (due  to  the  term  %R),  nonlinear  system  of  second-order  PDE,  with  the  flavor 

of  a  coupled  system  of  first-order  Hamilton-Jacobi  equations  involving  wave-phenomena, 
singularly  perturbed  by  a  second-order  term  (Laplacian)  but  in  a  nonlinear  manner. 


4,  Shearing  of  a  bar 


Consider  a  problem  where  only  the 


U, 


12 


(30) 


component  of  Ue  is  non-vanishing  and  all  fields  vary  only  in  the  x,  direction.  So,  x3  is  the  axis 
of  an  infinite  cylinder  say  of  square  cross  section.  Let  jq ,  x2  be  orthogonal  directions  in  the  cross 
section.  Then  the  only  non-vanishing  components  of  the  Nye  tensor  a  are 

^132^12,3  =  ^11  =’>  ^11  =  — ^12,3  • 

We  assume  that  the  only  nonvanishing  stress  components  are 


12 


=  t(Uu;2  +u’tj 


=  T. 


21 


which  vary  only  in  the  x3  direction  and  so  equilibrium  (f  /  =  0)  is  identically  satisfied 


regardless  of  the  x3  variation  of  Ue .  Assume  that  required  shear  tractions  can  be  mobilized  on 
the  lateral  surfaces  of  the  cylinder. 

For  simplicity,  we  now  assume  that  curlWR  is  symmetric  and,  like  the  stress,  that  ( curlWfj 
is  the  only  nonvanishing  component  (up  to  symmetry). 
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Now, 


— -  =  -curl  (axV) 

where  V  is  given  by  (27).  Consider  the  tenn 

f-=^Jjr^k  ;  f  =  XTa. 

The  only  non-vanishing  components  of  T jrark  is  given  by 

^2rarl  ~^2\a\\  ~  ^12^11 

so  that  the  only  non-vanishing  component  of  /  is 

fi  =  S22\^2\a\\=~''^2\aW  (31) 

Now 

(  OCX  f  )^  =  £jka,jfk  —£\jZarjf  =  £nzar\fz  +  £i23arlfz 

{aXf)r2=~arJ. 3  (32) 

(axf  )rx=ar2fv 
So 

[cur! (ax  f)\  =  s.p  (a x  f)^  =  £jsi  (a x /)rl>j  +  eb2  (a x  f)^ . 

Therefore  the  only  non-vanishing  component  of  curl  (ax  f )  is 

[curl (ax  f)]u  =  sm  (ax  f)n3  =-(ax  f\23  =  (anf3 ) ,3 .  (33) 

Of  course,  under  the  assumptions  we  have  made,  one  can  use  similar  arguments  to  deal  with  the 
term  curl  WR  in  the  expression  for  the  dislocation  velocity.  The  remaining  tenn  to  deal  with  is 

X  (curia)7  a.  Since  an  is  the  only  non-zero  component  of  the  dislocation  density,  the  only 


survivor  in 


Ma)  =e.  a 

J  ij  jmn  in, m 


is  ( curl  a)  =  ejMa] , ,  —  an3.  Therefore,  the  only  non-zero  survivor  in 

lx  (curlaf  a)  —  esjr  l (curl  af  }  air 


(X  (curl  af  a  j  =  e321  ((curl  af  j  an  =  —a113an .  (34) 

Noting  (31)-(34),  it  can  now  be  deduced  that  the  evolution  equation  for  the  only  non-zero 
component  of  the  dislocation  density  tensor  is 


an 


B(au) 


au  —  a 


11  (^12  "f  IXr)\2  ^  £an,f)  ’  • 

-  ,3 


We  now  assume  ul  2  to  be  the  only  non-zero  displacement  gradient  component  which  is  only  a 

function  of  x3  and  t.  Along  with  (30),  this  implies  that  U(’2  is  the  only  non-zero  plastic 
distortion  component  and  further  that 

an  ^i32^-212.3  ^2] 2,3  •  (36) 

Let  us  denote 
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<p\x,t)  Uf2  [x,t)  ;  x  =  x3  (37) 

Henceforth,  a  subscript  x  will  denote  partial  differentiation  w.r.t  x,  and  similarly  with  t .  Using 
(3 1),  (32)  2 ,  and  (34),  (29)  now  reduces  to 


<P, 


M2 
kl) 


B 


■12 


+  £<P*c+Zr12]> 


(38) 


using  the  fact  that  —div  grad x  =  curl  a  .  The  primary  reason  why  the  1-d  scalar  version  of  the 
Laplacian  in  %  can  be  written  in  terms  of  the  full  plastic  distortion  {(pxx)  is  that  in  this  simplified 
problem  div(axV)  =  0 . 

For  the  sake  of  simplifying  analysis  in  a  special  case,  we  would  like  to  have  (38)  solely  as  an 
equation  for  cp .  The  assumption  on  the  spatial  variation  of  uA  2 ,  along  with  a  fixed  boundary 

condition  on  the  base  of  the  bar  and  the  fact  that  ux 3  has  to  vanish,  implies  that  ul2=g  be  a 

function  of  time  alone,  so  that  the  tenn  T(x,t)  =  T(g(t)  —  <p(x,t))  in  (38)  does  not  pose  a 


problem.  However,  the  tenn  {Xr)v  's  problematic  as  it  cannot  be  written  as  a  function  of  (p .  For 

the  sake  of  simplicity  in  this  special  problem,  we  will  make  the  approximation  that  it  is  simply  a 
function  of  cp  denoted  by  P .  Thus  the  governing  equations  we  consider  are 


m 


m 


*(kl) 

M2 

kkl) 


r  {s  (0  -  <p)  +  £(Pxx  +  P{<P)\  hard  loading 
r  (x,  t)  +  £(pxx  +  P(<p)]  soft  loading. 


(39) 


Equation  (39)  governs  the  evolution  of  (generally  incompatible)  plastic  shear  strain  in  the  bar; 
spatial  gradients  of  the  profile  represent  the  dislocation  density  field  ((36)-(37))  with  individual 
discrete  dislocations  (more  precisely,  dislocation  walls)  represented  by  a  sharply  localized  tpx 
field  or,  alternatively,  a  sharp  transition  front  in  the  cp  field. 

Under  hard  loading  (displacement  control),  g(t)  represents  the  average  applied  engineering 
shear  strain  on  the  bar,  that,  under  the  controlling  assumptions,  can  only  vary  at  most  with  time. 
Under  soft  loading  (traction  control),  r[x,t)  represents  the  applied  shear  stress  on  the  top  surface 

of  the  bar  (of  course,  keeping  in  mind  that  constraints  would  have  to  be  in  place  that  can  provide 
the  required  reaction  tractions  on  the  lateral  surface  of  the  cylinder;  admittedly,  this  is  not  a 
practically  realistic  requirement). 

The  function  gif)  is  defined  as  follows:  by  assumption,  let  u12  =g0(x3,t)  where  g0  is  any 
function  of  the  arguments  shown.  Then,  u1(x1,x2,x3,t)  =  g0(x3,t)x2 +h(xt,x3,t),  where  li  is 
another  function  of  the  arguments  shown.  We  now  assume  that  the  base  of  the  bar  is  at  x2  =  0 
and  is  held  fixed,  i.e.  ux  (jc150,x3,?)  =  0  =>  h  =  0  .  The  applied  displacement  at  the  top  of  the 
cylinder,  x2=l,  is  given  as  g0  (x3,t)l ;  but  we  need,  by  assumption,  that  ux3  =  0.  Then 
g0  [x3,t)  =  g{t),  where  g  has  the  physical  meaning  mentioned  in  the  last  paragraph. 
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As  for  boundary  conditions,  we  either  consider  the  Cauchy  problem  on  the  domain  (—00,00), 
or  with  Nuemann  conditions  cpx  =  0  at  both  ends  of  a  finite  bar. 

We  note  that  this  problem  has  been  intentionally  set  up  to  be  simple  enough  so  that  force 
equilibrium  is  trivially  satisfied  and  hence  does  not  play  a  role  in  the  plastic  strain  evolution,  the 
latter  being  directly  controlled  by  the  applied  boundary  conditions.  In  fact,  stress  fields  arising 
from  any  solution  (39)  is  a  genuine  static  solutions  of  balance  of  linear  momentum  without  any 
approximation  on  inertia. 

4.1  P{(p)  =  0,  hard  loading;  exact  solutions  for  equilibria  in  a  special  case  and  formulation  of 
questions  of  stability 


We  deal  with  (39)  j  under  the  assumption  that  g  is  an  arbitrarily  fixed  constant,  i.e.  the  time 

dependence  of  the  loading  is  not  considered. 

The  equation  is 

<Pt  =  +  r(g-<p)],  (40) 

D 

where,  from  (27)  2  candidates  for  F  are 


F{a)  =  a 2  or 
F[a)  =  |a| 


(41) 


and  we  abuse  notation  and  refer  to  both  B,B*  as  B  . 


We  define  the  function  z  as  follows:  let  z  be  the  following  ‘periodic-cubic’  function: 


extended  periodically  beyond 


ju  >  0 . 


(42) 


Because  of  natural  physical  considerations,  we  need 

r(0)  =  0  ;  p  =  r'(0)  >  0  (43) 

i.e.  shear  stress  at  no  elastic  strain  should  vanish  and  the  elastic  modulus  at  zero  strain  should  be 
positive.  Now  define  z  in  (40)  as 

(44) 

Question  0:  Are  there  spatially  varying,  ‘dislocation  like’  equilibria  of  the  equation  (e.g. 
tanh(x)  type)  for  special  values  of  g  (or  a  family  of  spatially  inhomogeneous  equilibria 
parametrized  by  g )?  Of  course,  there  is  a  1  parameter  family  of  spatially  homogeneous 
equilibria. 

Inhomogeneous  equilibria:  It  can  be  checked  that  the  one-parameter  family  of  profiles 

<Peq{x>g)  =  g~  f  +  (45) 

2  2  (  V  4s  J 
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are  equilibria.  To  see  this,  define  ye'.=  g  —  (p  and  note  that  ye(x)  —  p/2  corresponding  to  the 
profiles  in  (45)  belong  to  the  range  [—  (s>/2,^/2] .  Next  note  that  the  periodic-cubic  defined  in 
(42)  agrees  with  the  actual  cubic  function  of  the  same  form  on  the  interval  [—  cp  /  2,  cp  jf\ .  Let  us 
call  this  ‘actual’  cubic,  i.e.  without  periodic  extension  as  in  (42),  c(-) .  Now  for  0  <  ye  (jc)  <  (p  , 
which  is  the  range  of  values  covered  by  the  elastic  strain  ye  of  the  profiles  in  (45),  a  calculation 
for  the  profiles  (45)  shows  that  £{(peq)  (x)  =  c[ye  (x)  —  #>/2) ,  for  each  x  in  (—00,00).  But  for 
0  <ye<(p,  c{ye  —  (pl2)  =  T(ye  —  (pl2)  =  —T(je)  by  construction  of  c  and  r  and  (44),  which 
implies  s{<P^)  (x)  +  r  (^e  (x))  =  0  for  all  x  in  (—00,00)  for  the  profiles  (45).  Of  course,  by 

design,  r  satisfies  the  physical  requirement  (43).  The  problem  definition  and  special  solution 
have  been  achieved  by  inspection:  we  consider  admissibility  of  equilibria  of  the  physically 
desired  form,  i.e.  A  +  B  tanh  (sx)  ,  ( A,B,s  constants),  and  pick  the  constants  and  design  the 
function  r  to  satisfy  physical  requirements  and  equilibrium. 

Now  consider  the  equilibrium  for  g  =  g„.  Consider  the  PDE  with  g(x,t)  =  gt  +  8,8  >  0  ,  where 
g, ,  8  are  constants,  and  initial  condition 

^(x,0)  =  cpeq  (x;g*) .  (46) 

Question  1:  The  question  is  whether  there  exist  travelling  wave  solutions  (with  non-zero 
velocity,  and  hopefully  “dislocation-like’)  for  equation  (40)  with  constant  g  =  g„+8  and  if  the 
initial  condition  (46)  is  attracted  to  such  a  travelling  wave  as  t  — >  00  or  to  other  equibiria  that 
exist,  in  particular  (peq  (x;g*  .  For  gt  =  0  if  such  a  travelling  wave  exists  for  some  Sp>  0  to 

which  (peq  (x;0)  is  attracted,  then  we  would  have  shown  the  existence  of  a  ‘Peierls  strain’  in  the 
model. 

Traveling  waves:  We  call  a  function  of  the  form  d(x,t):=  f(x  —  ct ),  c  a  constant,  a 
travelling  wave  solution  to  (40)  if  it  satisfies  it  with  c  ^  0 .  Here  /  is  a  function  of  a  single 
variable.  On  considering  the  governing  equation  of  a  traveling  wave  solution 

-f'c  =  H^p-{ef"  +  T(g,  +  <?-/)}  (47) 

n 

it  can  be  checked  that  (peq  (x  —  ct,g*^  is  not  a  travelling  wave  solution  to  (40)  with  g  =  g*  +S , 
for  8  ^  0,  |  <C  1 .  In  this  sense,  it  may  be  said  that  the  model  has  a  Peierls  barrier,  in  that  these 

‘dislocation’-like  equilibria  cannot  be  rigidly  translated  by  changing  the  load.  Question  1  above 
addresses  the  question  of  whether  there  exists  load  levels  corresponding  to  which  there  are 
profiles  that  move  rigidly  and  whether,  upon  additional  loading,  the  equilibrium  profiles 
displayed  above  in  (45)  can  be  defonned  into  ones  close  to  these  travelling  waves. 

Question  2:  (pt  (x, /)  =  0  is  an  equilibrium  solution.  Is  there  a  value  of  g  >  0  in  (40),  say  gN  , 

and  further  conditions  (e.g.  boundary  conditions)  for  which  <pt  becomes  ‘unstable’  and  some 
initial  condition  (pt  (x)  +  8(p(x)  (small  perturbation)  is  attracted  to  the  equilibria  (p<y  (x;  g  v )  ?  If 
so,  then  we  would  have  shown  the  possibility  of  nucleation  within  the  model  as  a  question  of 
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dynamic  instability  (cf.  Dayal  and  Bhattacharya  (2006)  in  the  context  of  peridynamic  theory). 
The  analogous  question  for  g  time-dependent  can  be  similarly  formulated. 


As  an  interesting  aside,  it  is  important  to  realize  that  a  time-parametrized  family  of  equilibria 
defined  by 


<Peq(x,t)  =  g{t)~j±j  tanhU— x 


where  g(-)  is  any  non-constant,  time-dependent  loading  program,  is  not  a  solution  to  the 
governing  equation,  for  the  same  time-dependent  loading,  given  by 

9,  =-^j^[£<Pxx  +  T(g(t)-<p)\-  (49) 

Such  are  the  curiosities  of  ‘quasi-static  solutions’  of  evolution  equations,  that  may  only  be 
justified  in  an  asymptotic  sense  corresponding  to  loading  histories  representing  short  bursts  of 
rapid  loading  interspersed  with  long  periods  of  loads  held  constant. 


4.2  P(<p)  ^  0;  equilibria  for  hard  and  soft  loading  in  a  special  case 


Recall  that  P(<p)  is  a  crude  simplification  for  a  term  that  is  in  reality  a  component  of  the  non¬ 
local,  curl  projection  of  R  =  s(\a\l)dG/d% .  At  any  rate,  this  needs  constitutive  specification. 
Let  us  make  the  assumption  that 

P{(p):=-2v{g-(p),  (50) 

where  g  is  a  material-specific  strain  value.  For  analytical  progress,  let  us  assume  the  stress 
response  in  this  case  to  be  given  by 

*{/)  =  ?{/)  (51) 

where  z  is  defined  by  (42). 

Equilibria,  hard  loading:  For  the  case  of  hard  loading,  the  governing  equation  for  the  plastic 
strain  evolution  becomes 

<Pt  =^&l[^  +  f(g-^)-2f(g-^)],  (52) 

D 

where  F  is  defined  in  (41).  As  before,  we  consider  the  simpler  problem  of  (52)  with  a  constant, 
but  arbitrary,  value  of  the  load. 

It  is  clear  that  spatially  homogeneous  profiles  are  equilibria,  regardless  of  the  value  of  the 
applied  strain  g  .  Spatially  inhomogeneous  equilibria  can  be  written  down  for  (52),  stated  for  the 
special  value  of  g  =  g  ,  as 

g±  — tanh  J—x  .  (53) 

2  V  4s 

It  is  easy  to  check  that  the  functions  in  (53)  cannot  serve  as  travelling  wave  profiles  for  (52) 
stated  for  g  =  g  +  8, 8  ^  0,  |<?|  -C 1 . 

Question  3:  A  natural  question  that  arises  is  the  following  one  related  to  nucleation:  for  (52) 
stated  for  g  =  g  or  for  g  time-dependent,  does  the  homogeneous  solution  (p[x,t)  =  0  lose 
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stability  so  that  a  minor  perturbation  on  (p(x,t )  =  0  as  initial  condition  drives  the  solution  of  (52) 
to  the  profiles  given  in  (53)? 

Equilibria,  soft  loading:  For  soft  loading,  we  consider  the  governing  equation 

<Pt  =^^\j{x)  +  e<pxx-2f(g-q>)'\,  (54) 

n 

where  r(x)  is  the  applied  traction  on  the  top  surface  of  the  bar. 

Again,  spatially  homogeneous  profiles  are  equilibria,  regardless  of  the  nature  of  r  .  For 
r  (jc)  =  0 ,  it  is  easy  to  check  that 


<  (*)  :=g±ytanh 


(55) 


are  equilibria  of  (54).  For  the  case  of  F(a)  =  |o| ,  it  is  interesting  that  qf  ( x  —  ct )  is  a  travelling 

wave  solution  to  (54)  stated  for  r  a  non-zero  constant.  Noting  that  (55)  represents  spatially 
monotone  profiles,  i.e.  the  sign  of  the  first  derivative  is  either  positive  or  negative  at  all  x ,  the 
constant  speed  of  the  travelling  wave  is  given  by 


c 


-sgn \{<P( 


eq 


r 

B 


(56) 


Thus,  this  case  corresponds  to  the  classical  conclusion  of  Peierls  (1940)  of  a  continuum  theory 
not  having  a  Peierls  barrier.  However,  as  is  clear  from  our  model  and  analysis,  translational 
invariance  of  the  theory  is  not  a  sufficient  condition  for  the  absence  of  a  Peierls  barrier  in  this 
dynamical  point  of  view. 

Dynamical  stability  questions  analogous  to  the  P[cp )  =  0  case  can  also  be  formulated  in  these 


cases  for  P(cp)  ^  0. 


5.  Implications  for  Meso/Macroscopic  Plasticity 


The  question  we  now  address  is  what  information  can  be  inferred  about  the  evolution  of  plastic 
distortion  averaged  over  space  and/or  time,  given  the  microscopic  evolution  equation  (29).  From 
(28)-(29)  we  note  that  the  plastic  distortion  evolution  may  be  written  in  the  form 


Up  = 


H\a\) 


axX(Sa), 


(57) 


where  S  is  an  appropriate  second-order  tensor.  Now, 

axX(Sa)  =  (aaT^ST  —  a S a . 


(58) 


We  assume  that  on  solving  the  microscopic  theory  and  probing  the  dislocation  density  field 
around  any  point  x ,  the  a  can  be  expressed  in  the  fonn 
a{x,t )  =  \a(x,t)\d(x,t)D(x,t) 


N(x,t) 

D(x,t) ^2  m1  (x,t),  d 

7=1 


1/|Z>|  if  D^O 
0  otherwise 


(59) 


where  d  is  an  0{\)  non-dimensional  numerical  factor,  N(x,i)  is  an  integer  that  represents  the 
number  of  dislocation  segments  involved  in  a  possible  junction  at  the  location  (jc,t) ,  m1  (x,t)  is 
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a  unit  vector  representing  the  Burgers  vector  of  the  Ith  dislocation  segment  that  may  be 
involved  in  the  junction,  and  similarly  l1  (x,t)  is  the  line  direction  of  the  same  segment.  Then 


aa 


=  a~  d  D  D,  8.  Si. 

ip  kp  jr  kj 


aSal  -  a' d  DiPkrSkj 


(60) 


S  =  T  +  [-£  div  grad  %  +  xR\  , 
and  defining  fourth-order  structure  tensor  fields 

V}^t)\rkj:=D'p^t)D^xB)5f 

the  evolution  of  (asymmetric)  plastic  distortion  takes  the  form 


2 

FTP-  a 

_ /-/  ^  ^  1  T _ p  arnr!  v  —1—  v  ^ 

-  ~  / 

a 

^  6/  y/L j  /Li  j .  yi  &  ciiv  grcici  %  i  j 

(61) 


(62) 


To  define  spatially  averaged  behavior,  we  resort  to  the  use  of  ‘filters.’  For  a  microscopic  field 
/  given  as  a  function  of  space,  we  define  the  meso/macroscopic  space  averaged  field  /  as 
follows: 

1 


/  w(x  —  x'\dx 

Jn(x)  v  > 


-J  w(x  — x'^f(x' ,t}dx' , 


(63) 


where  B  is  the  body.  In  the  above,  0[x)  is  a  bounded  region  within  the  body  around  the  point 
x  with  linear  dimension  of  the  order  of  the  spatial  resolution  of  the  macroscopic  model  we  seek. 
The  averaged  field  /  is  simply  a  weighted,  running  space-average  of  the  microscopic  field  / 
over  regions  whose  scale  is  determined  by  the  scale  of  spatial  resolution  of  the  averaged  model 
one  seeks.  The  weighting  function  w  is  non-dimensional,  assumed  to  be  smooth  in  the  variables 
x,x'  and,  for  fixed  x,  have  support  (i.e.  to  be  non-zero)  only  in  Q(x)  when  viewed  as  a 

function  of  (V) .  Applying  this  operator  to  (62)  (now  phrased  in  terms  x'  ),  we  obtain 
d  ^ 


~^-Up  =  Up  =  tn  d2  —  Z2J :  [T  —  e  div  grad x  +  Xr) 


under  the  assumption  that 
where 


UJ  ~  cu 


a 


57  = 


(64) 

(65) 

(66) 


The  assumption  (65)  is  valid  when  there  is  a  large  number  of  dislocation  segments  within  the 
averaging  volume.  It  is  also  assumed  that  the  average  magnitude  of  a  within  a  dislocation  core 
is  the  same  at  all  locations  where  a  dislocation  segment  is  present. 

The  physical  dimensions  of  B  is  given  by  =  .  Let  the  average  stress  magnitude 

in  the  core  of  a  single  dislocation  be  ~  cr* .  The  average  magnitude  of  the  dislocation  density  in 
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the  core  is  -bib1 .  According  to  Nabarro  (1987,  p.  506-507)  the  free-flight  velocity  from 
experiment  for  a  single  dislocation  is  vexp  =  \m/s.  -  \00m/s.  Thus,  the  magnitude  of  B  is 

~  a* / bvexp  and  the  magnitude  of  uj  is  ~  vexp  j bo  .  Let  C  be  the  time-resolution  of  observations 
and  define  5  =  t/t* .  Then  (64)  may  be  rewritten  as 

j-Up=  t* infer*  |  j-^d2  [z1  —  Z2] :  (T  —  sdiv grad X  +  Xr ) j  (67) 

where  /  is  some  non-dimensional  constant.  Let  us  define  the  nondimensional  viscosity,  v ,  by 

—  (68) 

v 

Assuming  the  magnitude  of  t*  and  b  to  be  l.v. ~  1 0  \y.  and  10  "bn.,  respectively,  and 
/  =  10~3  ~  1,  we  have 


10.1.10 

To 


^  1  ^  1.10  .1 

'  '  -10 


(69) 

10  1U  v  10  1U  ' 

Hence,  the  viscosity  v  «:  1 ,  and  (67)  takes  on  the  fonn  of  a  singularly  perturbed  system. 
Equation  (67)  implies  that  on  occasions  when  the  magnitude  of  the  spatially  averaged  plastic 
distortion  rate  on  the  t*  time-scale  is  not  too  large,  say  0(l/?*),  i.e.  at  most  ~1 /t* ,  the 
constraint 


Y  :=  d1  [Z1  —  Z2J :  [T  —  sdiv grad %  +  Xr)  ~  0 ,  (70) 

has  to  be  satisfied,  and  it  may  be  thought  of  as  a  tensorial  yield  criterion  governing  rate- 
insensitive  plastic  response  at  the  macroscopic  scale.  It  is  worthy  of  note  that  our  consideration 
does  not  state  that  the  constraint  (70)  has  to  hold  at  all  times;  indeed  it  allows  for  high-rate 
deformations  where  the  yield  condition  can  be  seriously  violated.  It  seems  to  be  certainly  within 
the  realm  of  possibility  that  Up  evolution  on  the  t*  time-scale  corresponding  to  variations  in 
applied  boundary  loads  shows  two  different  regimes;  one  very  slow  to  vanishing,  and  the  other 
an  order(s)  of  magnitude  more  rapid,  but  nevertheless  both  regimes  being  at  most  \/t*  .  Were 
this  to  be  the  case,  macroscopic  yielding  from  elastic  response  would  have  been  predicted, 
providing  a  fundamental  basis  for  qualitative  features  of  phenomenological,  macroscopic,  rate- 
independent  plasticity  theory.  Such  studies,  along  with  the  prediction  of  episodes  of  onset  and 
arrest  of  fully  rate-sensitive  response  with  significant  deviations  from  the  yield  set  would  be 
practically  useful  results4;  however,  they  clearly  require  a  deep  understanding  of  the  stability 
properties  of  the  yield  set  viewed  as  a  ‘limit  set’  of  trajectories  of  the  dynamics  (67)  with  all 
external  influences  held  constant.  For  a  glimpse  of  what  is  possible  in  situations  that  can  only  be 
considered  ‘simple’  in  the  context  of  our  general  theory  see,  e.g.  Carr  and  Pego  (1989)  and,  in 
the  context  of  ODE  theory,  Artstein  and  Vigodner  (1996),  Artstein  (2002)  (the  latter  has  low¬ 
dimensional  illustrative  examples  of  coupled  slow-fast  motions  with  non-trivial  limit  motions). 

Since  evolution  via  (67)  is  naturally  viewed  in  the  t*  time-scale,  the  condition  (70)  may  be 
viewed  as  valid  continuously  over  periods  of  the  order  of  t* ;  consequently,  (70)  may  also  be 
viewed  as  valid  when  averaged  over  a  period  of  t* .  Assuming  the  incompatible  slip  energy  G 


4  Preliminary  numerical  computations  indicate  that  both  yielding  type  response  and  intermittency  are  within  the 
predictive  scope  of  the  model,  at  least  in  the  1-d  case. 
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can  be  characterized  from  atomistic/quantum  studies,  characterization  of  the  time-averaged  form 
of  (70)  from  theoretical  homogenization  or  lattice-statics  calculations  or  simulations  of  the 
microscopic  theory  described  in  section  3  seems  to  be  a  worthy  and  practically  useful  challenge 
in  defining  macroscopic  plastic  response.  When  replacing  pointwise-field  values  at  (x',tr)  of  the 


individual  terms  appearing  in  (70)  by  their  space-time  averages  over  42  (x)  and  time  interval  t* 

around  the  instant  t  may  be  justified,  (70)  takes  the  instructive  form  (with  space-time  averages 
now  denoted  by  overhead  bars), 


Z‘-Z2 


T  -  £  div  grad  z  +  ZR 


=  0. 


(71) 


The  form  (71)  begins  to  suggest  possible  physical  origins  of  the  tenns  commonly  used  in 
phenomenological  plasticity  for  work-hardening,  the  “gradient  effects”  of  Aifantis  (1984,  1987)5, 
and  plastic  anisotropy  and  spin.  From  this  point  of  view,  it  is  interesting  to  note  the  nonlocality 
inherent  in  work-hardening  response.  Curiously,  if  one  were  for  the  moment  to  assume  the 
leading  structure  tensor  term  in  (71)  to  be  the  fourth-order  identity,  then  the  skew  symmetric 
part  of  the  equation  poses  an  interesting  integro-differential  constraint  on  the  skew  part  of  z  ■ 

We  note  here  that  the  association  of  high  gradients  of  slip  in  a  core  along  with  the  presence  of 
many  cores  as  the  primary  cause  for  rate-insensitive  macroscopic  response  is  different  and 
complementary  to  other  mechanisms  for  the  phenomena  based  on  pinning  due  to  postulated 
rugged  energy  landscapes  (Puglisi  and  Truskinovsky,  2005)  or  extrinsic  obstacles  (e.g. 
Bhattacharya  (1999),  in  the  context  of  phase  boundary  pinning).  It  is  also  interesting  to  note  that 
the  assumption  of  a  high  dislocation  content  within  the  averaging  volume  leads  to  wave-like 
response  completely  dropping  out  of  macroscopic  response  (when  the  average  plastic  distortion 
rate  is  small).  At  the  mesoscopic  scale,  this  key  assumption  of  a  uniformly  high  dislocation 
density  in  space  cannot  be  made;  then,  the  transport  terms  remain,  but  the  microscopic  theory  is 
practically  untenable.  What  the  form  of  the  governing  equations  should  be  at  this  scale  is  a  key 
challenge  for  rigorous  mechanics  and  mathematics.  Some  promising  heuristic  advances  coupled 
with  numerical  work  have  been  made  recently  (Acharya  and  Roy,  2006;  Roy  and  Acharya,  2006; 
Taupin  et  al.,  2007;  Acharya  et  al.  2008;  Taupin  et  ah,  2008;  Fressengeas  et  ah,  2009). 


6.  Comparison  with  the  structure  of  phase-field  models  and  level  set  evolution 

We  end  with  a  couple  of  remarks  on  the  structure  of  our  theory.  While  the  manner  in  which  the 
constitutive  structure  (27)  is  defined,  along  with  the  assumption  (13)  and  the  use  of  a 
conservation  law,  i.e.  (10),  may  be  reminiscent  of  (time-dependent)  Ginzburg-Landau,  local- 
energy-minima  seeking  phase-field  type  models  or  Cahn-Hilliard  models,  the  dynamics  implied 
by  our  theory  is  different  from  that  implied  by  these  fonnalisms.  Roughly  speaking,  this  is 
because  our  situation  corresponds  to  having  a  conservation  law  for  the  gradient  of  the  order 
parameter  (plastic  distortion)  based  in  the  kinematics  of  convective  transport,  along  with  a 
‘surface  energy’  term  -  alternatively,  if  one  considers  the  dislocation  density  as  the  order 
parameter,  then  there  is  a  conservation  law  but  no  surface  energy.  To  see  one  of  these 
comparisons  in  the  simplest  possible  context,  consider  an  energy  of  the  form 


5  Although  beyond  the  Laplacian  and  the  sign  in  front  of  it,  the  details  are  quite  different.  Note  that,  roughly,  Z 
corresponds  to  —  Up  . 


17 


(72) 


1  2 

l//  =  lj/{g-(p)  +  -{(Px)  , 

where  g  is  a  constant.  Then  the  corresponding  Ginzburg-Landau  model  would  be 

cp  =  ~r  r  =  mobility 

&P  (73) 

^<P  =  r(<Pxx  +  T{g-(p))’ 

and  it  is  instructive  to  compare  this  with  (40).  Of  course  the  all  important  difference  is  the 
F  ( cp  )  term  in  our  dynamical  model,  one  of  the  major  effects  of  which  is  in  providing  the 

singular  perturbation  structure  in  the  averaged  macroscopic  model  of  plasticity  mentioned  in 
Section  5.  This  comment  is  also  the  only  relevant  one  we  are  able  to  make  in  a  comparison  to  the 
static,  energy  minimizing,  phase  field  framework  of  Koslowski  et  al.  (2002),  where  it  is  not  clear 
exactly  what  physical  criteria  is  used  to  pick  amongst  multiple  local  minima.  Also,  it  is  well- 
known  that  the  Cahn-Hilliard  formalism  yields  fourth-order  equations  whereas  ours  is  simply 
second  order. 

Finally,  we  mention  that  our  equation  (39)  (with  B  =  B\tpx\)  in  the  scalar  case  may  be 

interpreted  as  the  evolution  of  a  level-set  with  velocity  given  by  the  Ginzburg-Landau  driving 
force.  Thus  it  is  natural  to  expect  transfer  of  qualitative  properties  like  pinning,  or  lack  thereof, 
between  our  model  and  level-set  propagation  of  this  specific  type. 
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